Indistinguishable Photons from a Single Molecule 
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Abstract 

We report the results of coincidence counting experiments at the output of a Michelson inter- 
ferometer using the zero-phonon-line emission of a single molecule at 1.4 K. Under continuous 
wave excitation, we observe the absence of coincidence counts as an indication of two-photon in- 
terference. This corresponds to the observation of Hong-Ou-Mandel correlations and proves the 
suitability of the zero-phonon-line emission of single molecules for applications in linear optics 
quantum computation. 
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Two otherwise indistinguishable photons arriving at a 50/50 beam splitter from different 
input channels, will interact and leave through the same output channel jO]. This phe- 
nomenon of two-photon interference is due to the bosonic nature of photons. Its observation 
is an important demonstration of the quantum theory of light. Apart from its fundamental 
aspect, two-photon interference also lies at the heart of a broad range of applications in 
quantum information science. Photon-photon interactions mediated by two-photon inter- 
ference allow for quantum gate operations employing photonic qubits. In such a scheme, 
linear optical elements, ideal photon counters, and true single photon sources were proven 
to be the only elements necessary to perform quantum computation [2] . The requisite single 
photon sources for these applications should provide indistinguishable (i.e. exhibiting ideal 
two-photon interference) single photons on demand. In practice, the construction of such 
single photon sources is one of the major tasks for the realization of linear optics quantum 
computation schemes. To date, experimental demonstrations relying on two-photon interfer- 
ence mostly used photons generated in a parametric down conversion process jjj. However, 
with this photon source, the probability of simultaneous two-photon emission is governed 
by Poissonian statistics. Therefore, only low photon fluxes can be used. This limitation has 
been one important motivation for the recent development of a variet y of true single photon 
sources based on different sing le two-level emitters fl Q, Q, S S HQ • These single pho- 
ton sources can emit single photons with a sub-Poissonian photon statistics provided that 
a large collection efficiency is realized. Only very recently, the coherence properties of the 
emitted single photons have been investigated. Indistinguishability has been demonstrated 
for photons emitted by single quantum dots Q,[l2, 1A\ and single trapped atoms [l^ . 



Single molecules are an attractive alternative to both, single quantum dots and single 
trapped atoms. In self-assembled InAs quantum dots, the achievable coherence length is 
limited to ~ 750 ps Q, Il2|- By contrast, dephasing in single molecules is mainly due to 
phonon coupling which is suppressed at cryogenic temperatures. Nearly transform limited 
emission from single molecules with coherence lengths close to 4.9 ns is therefore readily 



observed using vibronic excitation at 1.4 K 



In comparison to single atoms, single 



molecules offer the advantage of much longer observation times. Experiments with one single 
molecule over more than 20 days have been reported [19]. Single molecules have previously 
been employed as a source for triggered single photons using rapid adiabatic passage 5( 



and pulsed optical pumping |4J, Together with the large observed coherence times, these 
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experiments prove the potential of single molecules as sources for indistinguishable single 
photons. In this letter, we report Hong-Ou-Mandel experiments using the zero-phonon- 
line (ZPL) emission from a single terrylenediimide (TDI) molecule. As a result of these 
experiments, two-photon interference is demonstrated. 

Recently reported experiments on two-photon interference using single quantum dots 
111] . Q and single trapped atoms |3] all employed pulsed excitation schemes. In these 
cases, a Michelson interferometer with a path length difference that is exactly equal to 
the pulse separation is employed. This guarantees that consecutively emitted photons will 
simultaneously arrive at the two inputs of the beam splitter with a non- vanishing probability. 
Hence, depending on the coherence properties of the photons, two-photon interference is 
observed. In contrast, we have employed continuous wave (cw) excitation in our experiment. 
Under cw excitation, the path length difference between the two arms of the Michelson 
interferometer is not matched to an exact value. Rather, a path length difference which is 
larger than half the coherence length of the emitted photons is selected. This is sufficient to 
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ensure that two independent photons with negligible magnitude squared overlap integrals 
will simultaneously arrive at two inputs of the beam splitter with a certain probability. The 
signature of two-photon interference is then revealed as the lack of coincidence counts at the 
beam splitter output. 

The Michelson interferometer used in our experiments is depicted in Figure 1. We have 
used a rotatable A/2 plate in the longer arm of the interferometer in order to achieve parallel 
or orthogonal polarizations in channels 1 and 2. Considering parallel polarizations, a coin- 
cidence counting experiment between channels 3 and 4 reveals the normalized second-order 
coherence function: 

(2) = (al(t)al(t + r)o 4 (t + r)o 3 (t)) 
^ 34l ' Tj (ai(t)a 3 (t))(ai(t + T)a i (t + T)) , 

r represents the delay time between channels 3 and 4. In general, the expression of #34 (£, t) 
in terms of a, the photon annihilation operator of the molecular emission at the entrance 
of the interferometer, will be composed of 8 correlation functions and the corresponding 
complex conjugates. For this general case, the single photon wave functions in channels 1 
and 2, described as \ip to ) oc dte~ rspo " ( '~' o)/2+ ^ ( *M(t)|0), are not orthogonal. t denotes 
the initial time of the photon wave packet and <p(t) is a random function describing the 



pure dephasing process j20|,|2lj. The magnitude square of the inner product of delayed and 



undelayed single photon wavef unctions, | l^ — At) | 2 , can be shown to be proportional to 
e _27A * where 7 = 1/T 2 = T spon /2 + 7 pure is the total dephasing rate of the ZPL including 
spontaneous emission (r spo „) and dephasing due to other sources impure)- 

For our case of interest, in the limit of large At (At 3> 1/(27)) an d small r (r <C At), 
Eq. ^ is largely simplified. In this limit, the two-photon wave function in channels 1 and 
2 can be described as an outer product of two independent single photon wave functions 
iVOlV'to-At) with I {^to l^fo-Ai) | 2 = 0. Thus, the delayed and undelayed operators will satisfy 
the commutation relationship: [a(t — At + Ti),aJ(t + r 2 )] = 0, where 77 and r 2 represent 
small delays {ti,t 2 < 7) in the coincidence experiment. In the case of parallel polarizations 
at the exit of the interferometer, gf} (t,r) is then expressed as [22! I^: 



9% (r) = £ (V 2) (r) + 1 j - "" 4 Vr: 4" 9 W W ' (2) 
11 2 \ > cos 4 1> + sin a 

In this equation gW( T ) = and g^{r) = ^fe^^ correspond to the 

normalized first-order and second-order coherence functions of the ZPL emission of a single 

molecule respectively. Note that we drop the dependence of the coherence functions on t due 

to cw excitation conditions. Transmission and reflection coefficients in the beam-splitter are 

noted as cos 2 9 and sin 2 6 respectively. For mutually orthogonal polarizations, photons at 

both input channels are completely distinguishable. For this case, the result of a coincidence 

experiment between channels 3 and 4, reveals: 

&(r) = U^i^ + l). (3) 



1 //■«*, , \ sin 2 9 cos 2 
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Vibronic excitation as used here relies on the excitation of a fast relaxing (relaxation 
time ~1-10 ps) high energy vibrational level. In such a three-level incoherent excitation, 
the first and the second order coherence functions are given by analytical expressions as 

gW(r) = e~ 7T and g( 2 \r) = 1 — e ( w >- r V») T W P corresponds to an effective pumping 

(2) 

rate. Using these equations, an examplary solution of g^ 4 (r) is depicted for parallel (solid) 
and orthogonal (dashed) polarizations. As Eqs. El H3 assume small delay times, Fig. 2 is 
plotted for delay times between — 0.5/T spon and 0.5/r spon . In the solid curve, the signature 
of two-photon interference is the absence of coincidence events around zero delay time. 

The vibronic excitation scheme is depicted in Fig. 3a Q, l^. In short, we used a 
narrow band cw laser excitation tunable around 605 nm (Coherent 899) to excite a vibronic 
transition of TDI embedded in the Shpolsk'ii matrix hexadecane. The samples were prepared 



by adding hexadecane to a solution of TDI in CHCI3. CHCI3 and O2 were removed by several 
freeze and thaw cycles using liquid N 2 and intermediate evacuation steps. After this, the 
TDI/hexadecane solution was saturated with Ar. A drop of this solution was then quickly 
inserted into the precooled cryostat. Finally, the sample was cooled down to 1.4 K in a 
liquid He bath. An aspheric lens (NA=0.55) was used to focus the excitation and collect the 
emission light. The emission from the ZPL of a single TDI molecule was then filtered with 
a narrow band (FWHM = 1 nm) interference filter before being focused onto a confocal 
pinhole (75 \xm diameter) placed inside a 4x telescope. In order to select single molecules 
from different spatial positions in the sample, a piezo scanner (Attocube ANP100) was used. 
Once a molecule was selected, its emission could be observed over several days. Fluorescence 
was either directed to a spectrograph (Jobin-Yvon HR460 with an Acton Research CCD 
camera, spectral resolution 30 GHz) or a homebuilt Michelson interferometer (Fig. 1). 
The path length difference between the two arms of the interferometer was equivalent to a 
temporal delay of At = 4.6 ns. One of its arms was equipped with a rotatable A/2-plate. 
The output signals at the beam splitter were detected using avalanche photodiodes (APD, 
Perkin Elmer) in a Hanbury-Brown and Twiss configuration. The APD outputs were sent 
to a time to amplitude converter (TAC) as start and stop pulses. After adding an electronic 
delay on the stop pulse, the output of the TAC itself was stored in a multi channel analyzer 
(MCA). The resulting MCA histogram corresponds to the second-order correlation function 
#34 (t) which is normalized using the function's value at large delay times. A temporal 
resolution of 420 ps was achieved. 

A high resolution emission spectrum of the single TDI molecule is shown in Fig. 3b. A 
resolution limited band accompanied by the broad phonon sideband and vibronic emission 
lines is seen. The intensity ratio between the ZPL and all of its sidebands is governed by 
the Franck-Condon and Debye- Waller factors. In the emission spectrum in Fig. 3b, 40 % of 
the single molecule emission intensity is orginating from the ZPL. This ratio can be further 
improved by placing the molecule inside a cavity We investigated this molecule for 

five days including several cooling cycles between 1.4 K and 100 K. All data presented 
in this work have been recorded from this molecule. Occasionally, spectral jumps of the 
molecule were observed during the experiment (Fig. 3c). Excitation of the molecule at the 
new spectral position always allowed to bring back the absorption to the initial position at 
14988 cm" 1 within a few seconds. A similar switching behavior was reported previously [19J. 
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The digital spectral jumps are a clear signature of single molecule observation. In these rare 
events, the APD counts were observed to drop to < 5% of the original level. Hence a signal 
to background ratio > 95% was achieved in the experiments. 

The result of coincidence measurements in which the ZPL emission of the single TDI 
molecule was sent into the interferometer and recombined at the beam splitter without any 
polarization changes is depicted in Fig. 4a. In this normalized correlation function, at the 
zero delay position, the value of g( 2 \0) becomes 0.4. By contrast, Fig. 4b depicts the result 
of the coincidence detection experiment when rotating the polarization of the fluorescence 
in one interferometer arm by 90°. Clearly, the resulting correlation curve exhibits higher 
g( 2 '(r) values around zero delay time. The normalized difference between the two curves is 
shown in Fig. 4c. In this curve, the peak around zero delay time constitutes the proof for 
the observation of Hong-Ou-Mandel correlations. Figures 4d and 4e show the normalized 
difference between the data of Fig. 4a and another measurement with parallel polarizations, 
as well as for the data of Fig. 4b and another measurement with orthogonal polarization 
on the same molecule. In contrast to Fig. 4c, no signal above the noise level is seen at zero 
delay in both of these cases. 

Ideally, if the photons at the input of the beam splitter were completely indistinguishable, 
the value of g^(0) in Fig. 4a would be 0. Instead we observe g^ 2 \0) ~ 0.4 which corresponds 

9 < 2 ) (0)-g( 2) (0) 

to a coincidence reduction factor, V(0) = — - — > of V(0) = 0.24. The coincidence 

reduction factor can be lowered for several reasons. Imperfections in the mode matching of 
our interferometer are one source of contrast reduction. Because of its higher throughput, a 
conventional Michelson interferometer was preferred over a setup using a single mode glass 
fibre in our experiment. With a reduced interferometer path length difference of 1 ns, we 
observed a fringe contrast of 0.7 for single molecule emission. However, due to the loss of 
single photon interference, we were not able to determine the fringe contrast for the 4.6 ns 
path length difference used in the correlation measurements, but the value of 0.7 is a upper 
boundary for this case. In addition to interferometer imperfections, spectral dynamics below 
the 30 GHz resolution of our monochromator could deteriorate the contrast. If these jumps 
occur on the ns time scale of the excited state lifetime, they will give rise to dephasing and 
line broadening. While we did not determine the ZPL emission linewidth 7 of this molecule 
with a spectroscopic measurement, we showed previously that stable 65 MHz emission is 
readily achieved in this sample system 3]. Considering a typical spontaneous emission 



6 



rate T 

spon ~ 1/3.4 ns ^[18]) & fif to the correlation measurements in Figs. 4a and 4b yields 
Ipure ~ 1/5 ns -1 . Hence a photon coherence time of ~ 3 is derived. For this 
3 ns coherence time determined from our experiments and 4.6 ns path length difference, 
the magnitude squared overlap of the single photon wavefunctions in channels 1 and 2, 
KV'iolV'to-At)! 2 = e _27At , reveals 0.05 instead of 0. Thus, although our experiments clearly 
demonstrate the observation of Hong-Ou-Mandel correlations, Eqs. El and H3 are only valid 
to a good approximation due to the assumptions At ^> 1/(27) an d t <C At. A complete 
solution of Eq. ^ is necessary for a full theoretical explanation. It should be noted that 
a complete analysis will also explain the general appearance of the correlation functions 
depicted in Fig. 4 for large delay times (r > l/r spo „). Such an analysis should reveal the 
sidelobes apparent in Figs. 4a and 4b at r = ±At. 

In conclusion, we have demonstrated that the ZPL emission of single molecules can be 
used to produce indistin guis hable photons. This complements earlier experiments using 
single quantum dots jl^l ll^ and single trapped atoms [l^. Single molecules are an 
attractive alternative to these systems since they combine long coherence times with long 
observation times. We observed a coherence time of ~ 3 ns over the course of five days. 
The level of observed interference can be improved by coupling the single molecule emission 
into a glass fibre. This, however, necessitates an improvement of the collection efficiency. 
Experiments using specially designed low temperature optics with a high numerical aperture 
are currently pursued. Finally, the experiments presented here can easily be extended to 
include a pulsed excitation scheme. This is not necessary for the demonstration of two- 
photon interference but will be important for practical applications. 
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FIG. 1: Michelson interferometer used for two-photon interference measurements. 



FIGURE 1, Kiraz et al. 
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FIG. 2: Solution of <7 34 (r) for large At (At 3> I/7) and small r considering parallel and orthogonal 
polarizations. Parameter values are Wp = 2.5T spon , 7 pure = 3T spon , and = 7r/4. 



FIGURE 2, Kiraz et al. 
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FIG. 3: (a) Vibronic excitation scheme, (b) High resolution emission spectrum of the TDI molecule. 
Two different magnifications are plotted for clarity, excitation intensity is 4 mW. (c) Consecutive 
high resolution emission spectra under an excitation intensity of 1 mW. 



FIGURE 3, Kiraz et al. 
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FIG. 4: Measured normalized second-order photon correlation functions for parallel (a) and orthog- 

(2) (2) 

onal (b) polarizations, (c) Difference between the correlation functions g ± and normalized 
to g^ . Normalized difference between two measurements at parallel (d) and orthogonal (e) po- 
larizations. Excitation intensity is 1 mW in a-e. Binning of the data is 2 in a and b, and 7 in 
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FIGURE 4, Kiraz et al. 
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